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Abstract 

It has recently been argued that there might exist a four-parameter analytic solu- 
tion to the Bianchi IX cosmological model, which would extend the three-parameter 
solution of Belinskii et al. to one more arbitrary constant. We perform the pertur- 
bative Painleve test in the proper time variable, and confirm the possible existence 
of such an extension. 
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1 Introduction 



The Bianchi IX cosmological model results from the assumption of a universe with a 
homogeneous anisotropic space, and an important question concerns the singularities in 
the complex t-plane of the components (a 2 ,b 2 ,c 2 ) of the metric tensor [l], ^j, gfl. This 
system is made of three coupled second order nonlinear ordinary differential equations : 

2 a ( i d(lna)\ 4 2 2 2 

2a abc— abc = a — [b — c ) and cyclically, (1) 

at \ at J 

with the first integral 



a 2 



^ 2 d(b 2 ) d(c 2 ) + ^djc 2 ) dfy) + c2 d(a 2 ) d(6 2 



+ a 4 + 6 4 + c 4 - 2 (6 2 c 2 + c 2 a 2 + a 2 ^ 



dt dt dt dt dt dt 

(2) 

Here a 2 = is +1 or —1 according as the metric is asymptotically Euclidian or Minkowskian. 
The study of these singularities is easier in the so-called logarithmic time variable r, 
related to the proper time variable t through the transformation : 

dt = abc dr. (3) 

The main result of the Painleve test [Q— |J7| is a probable chaotic behaviour of the system. 
Moreover, one can prove || the inexistence of rational first integrals other than (0). As a 
by-product, this test can also be used to detect all possible particular analytic solutions, 
in addition to the three particular solutions which were already known. 

The first exact solution results from the existence of the Euler subsystem (1750) : 

ex — - — = lu 2 ^3 and cyclically (4) 
dr 

a ^° = a ana cyclically (5) 

with uji = be and cyclically, reducible to Weierstrass' elliptic equation in the variable 
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The three-parameter solution of system (|J), single- valued in r, was found by Abel and 
Jacobi, and later in the study of the Bianchi IX system by Belinskii et al. H|. 



A second subsystem exists, the so-called Darboux system [10 



a— — ^ = ^2^3 — ^1^2 — ^1^3 and cyclically (8) 
dr 

d(bc) 

a — - — = a — b — c and cyclically (9) 
at 
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which has been integrated by Halphen [11, 12]. This three-parameter solution, single- 



valued in r was rediscovered in the context of Bianchi IX, by Gibbons et al. ||13||. This 



system is equivalent to the third order Chazy equation of class III [[14] 



with y = — 2(o>i + uj 2 + u 3 )/a. 

Last, a third particular solution was found by Taub [|1|, when the metric tensor is 
axially symmetric (e.g. b 2 = c 2 ), leading to a four-parameter trigonometric solution. 

As to the Painleve analysis, it identifies three possible local behaviours of the general 
solution in the r variable 

1. a simple pole for a 2 and a simple zero for b 2 and c 2 ; 

2. a simple pole for a 2 , b 2 and c 2 ; 

3. a regular behaviour for a 2 and a double pole for b 2 and c 2 . 

Finally, the only possible new solution isolated by the Painleve analysis is || an exten- 
sion of Belinskii's elliptic solution to four parameters, whose closed form is yet unknown 
and which we denote for convenience BGPP4. The finding of this closed form has a strong 
interest in cosmology, since it would ipso facto yield a six-parameter single- valued solution 



to the Brans-Dicke Bianchi IX cosmo logical model |Tg, [17 1 



The aim of this paper is to further examine this four-parameter solution. The con- 
nection between the local behaviours in t and those in r is established in section |2[. We 
perform in section [^ the perturbative Painleve test Jl8|, 19] on system (|l|) in the proper 



time variable t, in order to see whether or not the existence of the four-parameter solution 
can be denied. 



2 The local study 

The successive steps are the following. First we determine all possible leading behaviours 
of the Laurent series for the dependent variables a(t),b(t) and c(t). This is done by 
substituting respectively a (t — t ) pi , bo(t — t ) P2 , co(t — to) P3 into the system (|l|), and by 
checking all the different cases for the constants p±, P3, under the condition that the 
constants a^, b , cq are all different from zero. 

Then, for each leading behaviour, we have to look for the so-called indices or resonances. 
These indices can be found by inserting : 

a{t) = a x Pl (l + ea rX r ); b(t) = b oX P2 (l + eb rX r y, c(t) = c oX p *(l + ec rX r )- 

with X = t — to, into the system (|l|), and by taking only those terms which are linear in 
e. This leads to a coupled system of three linear equations in a r , b r , c r , whose solution is 
unique iff the determinant of the coefficients differs from zero. The indices are the zeroes 
r of this determinant. 
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If all the leading exponents are positive, it is recommended to invert at least one of 
the dependent variables, for example a(t) a(t) _1 in order to calculate the indices. For 
a more detailed explanation we refer to f20fl . 

If we now apply the scheme explained above, we obtain only three possibilities, listed 
in Table |I[ 

The arbitrary values of the two opposite indices of the third family (III) immediately 
imply an infinite amount of logarithmic branching, hence probably chaos. 

The last row represents the local transformations, after integration of relation @ in 
which we introduced the leading behaviours for a(t), b(t) and c(t). 

Table 1: Families of singularities of a, b and c in the variable t. Notation is e 2 = e 2 = 1. 





(I) 


(II) 


(HI) 


(Pi, Pi, Pa) 


(-1 I i) 

V 3' 3' 3^ 


(1, 1, 1) 


(1, 0, 0) 


(oo, b , c ) 


(2eab c /3, b , c ) 


(e x , e 2 , £ 3 )/(2cr) 


(a , 6 , eb ) 


r 


-10 0--- 


-4,-4,-4,-1,2,2 


-1, 0, 0, 0, -2/(aa ), 2/(aa ) 


T - T 


(t - to) 2/3 


(t-to)- 2 


ln(t - t ) 



Analogously, Table || gathers the results of who performed the same kind of 

analysis in the r- variable and for the dependent variables A = a 2 , B = b 2 , C = c 2 . 

The comparison of the two tables immediately shows the specificity of family (III). 
Indeed, the local transformation between the two times is different for family (III). 

For families (I) and (II), the local transformation yields the correspondence between 
the exponents pj and the indices r. On the contrary, for family (III), the two local laws 
t — t = (t — to) -1 an d t — T = ln(t — to) ar e incompatible, and the deep reason for that 
is the non-Fuchsian nature of family (III) in the logarithmic time r. 
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Table 2: Families of singularities of A = a 2 , B = b 2 and C = c 2 in the variable r. 





(I) 


(II) 


(III) 


2(pi, Pa, p 3 ) 


(-1,1,1) 


(-1,-1,-1) 


(0, -2, -2) 


(A), -Bo, C ) 


(<T, So, Co) 


(cr, cr, a) 


(A), a 2 Mo, a 2 Mo) 


r 


-1,0,0,1,1,2 


-1,-1-1,2,2,2 


-1,0,0,2 


T - T 


(< - to) 2/3 


(t-t y 2 


(t-to)- 1 



3 Recovering the exact solutions 



Let us now perform the perturbative Painleve analysis [|18 , p79|| in order, mainly, to test 



the existence of the assumed four-parameter extension to the BGPP solution. 

Each of the three families is separately a local representation of the general solution, 
so each family must describe every closed form particular solution. 

For the family (I), we have checked the absence of any movable logarithm, so none of 
the known closed form particular solutions in r is excluded. 

For the family (III), the matrix P(r) of the usual recurrence relation 



Vr P(r) \b r +Q r = (11) 



P(r) = 2bl 




is 

/ a b a 2 r(r + 1) a^a 2 r a\ea 2 r 

a 2 o 2 r 2 - 2 2e j . (12) 

V 2 e(a 2 a 2 r 2 - 2) 

So the number of movable logarithms is : 2 (from the two irrational indices), plus 1 (from 
the difference between the multiplicity 3 and dim(Ker(P(r))) = 2 at index r = 0). When 
one returns to the r variable, the local transformation r — r = ln(t — t Q ) suppresses one 
of these three logarithms, leaving open the possibility of a four-parameter solution, single 
valued in the r variable. 

Let us now process the family (II), by expanding a, b, c in the double series 

+oo +oo +oo 

a = E £n E afk j+p = £ n a (n \ and cyclically, (13) 

n=0 j=np n=0 

where p stands for the lowest index p = —4. 

Due to the invariance under permutation, it is convenient to introduce the seven 
new variables (X, Y, Z, U, V, W, T) 
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b -l = T2 U+ 7S V + Ts W ( 15 ) 

a W 1 + b W 1 + c^ = T (16) 

where the six arbitrary variables are : (a-\,b_\,c_\) : two out of (c4 , 63 > c 2°^); an d one 
out of (a_], b_{, c^j). The necessary conditions for the absence of movable branching are 
the following. 

Order zero : n = 

We obtain at (n,j) = (0,2) the single condition Z = 0, i.e. + &2 + <4 = 0. Parity 
implies a§+i = 6§+i = c 2?+i = with j = 0, 1, 

Higher orders : n > 1 



At first order one must distinguish three cases, and examine them separately at higher 
orders. 

Case 1 : X = Y = 

Because a^\, b_\, cl\ and for instance ol\ are arbitrary, one can set for n > 1 
a^ 1 ], 6^4, cl} and equal to zero; the same can been done for and c% with 
n > 0. Each series a^ n ^ then terminates 

Vn : a <"> = t £ a^ 4n t 3j - 4n with a« = &3 = = i 

j=0 3 

and we have checked the absence of any logarithm up to perturbation order n — 8. 

This results in a four-parameter Laurent series depending on U, V, W, T. The three- 
parameter solution of Belinskii et al. corresponds to W = 0, as can be checked by substi- 
tuting the series (p~3|) into the Euler system (|5|). 

Case 2 : (W 2 - U 2 )U = 0, (3Y 2 - X 2 )X = 0. 

No additional no-log conditions are encountered at higher orders n, and this case 
isolates the three Taub reductions : b = c ^ a and cyclically. 

Case 3 : U = V = 0. 

At second order we find the only additional condition W = 0, which implies : 
af ] = bf> = cf ] = with j < 0, n = 1, 2, ... except for = 6^1 = c { l\ = T/3. 

This represents a three-parameter solution depending on X, Y, T, which is precisely the 
Halphen solution to the Darboux system @, as can be checked by direct substitution. 

The above results are summarized in Table [| together with their interpretation. 
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Table 3: Results of the perturbative Painleve test. 



order (n) 




n=0 


Z = 


n = 1 


X = Y = 


(3V 2 - U 2 )U = 
(3Y 2 -X 2 )X = 


U = V = 


n = 2 


no conditions 


(b - c)(c - a)(a - b) = 


= 


arbitrary parameters 


U, V, W, T 


X, U, W, T 


X,F,T 


identification 


BGPP4? 


Taub 


Halphen 



4 Conclusion 

We have applied the pertubative Painleve analysis to the Bianchi IX cosmological model 
in the proper time variable and connected the results to the already existing ones in the 
logarithmic time variable. The main result is the absence, checked to a high perturbation 
order, of any obstacle to the existence of a four-parameter extension to the solution of 
Belinskii et al, which would be single valued in the logarithmic time. 

However the challenge still remains to obtain the closed-form expression of this four- 
parameter extension. 
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